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Abstract 

We study the Kitaev-Ising model, where fen'omagnetic Ising interactions are added to the Ki- 
taev model on a lattice. This model has two phases which are characterized by topological and 
ferromagnetic order. Transitions between these two kinds of order are then studied on a quasi-one 
dimensional system, a ladder, and on a two dimensional periodic lattice, a torus. By exactly map- 
ping the quasi-one dimensional case to an anisotropic XY chain we show that the transition occurs 
at zero A where A is the strength of the ferromagnetic coupling. In the two dimensional case the 
model is mapped to a 2D Ising model in transverse field, where it shows a transition at finite value 
of A. A mean field treatment reveals the qualitative character of the transition and an approximate 
value for the transition point. Furthermore with perturbative calculation, we show that expectation 
value of Wilson loops behave as expected in the topological and fen^omagnetic phases. 
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1 Introduction 



In the study of many body systems, specially those inspked by quantum computation, a new paradigm 
is emerging, which embodies concepts such as topological order, topological phase, and topological 
phase transition. Contrary to the traditional Landau paradigm, topological phases are not character- 
ized by local order parameters and topological phase transitions are not accompanied by spontaneous 
symmetry breaking. In addition to their interest in condensed matter physicsfl], i.e. in fractional 
quantum Hall liquids |2i| and quantum spin liquids 13J, lattice models exhibiting topological order are 
of immense interest in the field of quantum computation and information ||4l, jS) due to their robust- 
ness against decoherence. The simplest such lattice models is the Kitaev model H, although other 
models like color codes have also been introduced and extensively studied 161 1?]. The ground state of 
the Kitaev model on a surface of genus g, exhibits a 4g-fold degeneracy which is directly related to the 
topology of the surface. Different ground states look exactly the same if probed by expectation values 
of local observable and are only distinguished if probed by global string-like operators going around 
non-trivial homology cycles of the surface. One can thus use these ground states to encode 2g qubits 
which are robust against errors and decoherence. Moreover one can do topological computations on 
these qubit states if one uses braiding and fusion of anyonic excitations of these models. 

It is then natural to ask how much this topological order in the original Kitaev model or its gen- 
eralizations to Zjv group or the topological color codes |i6j iTj are resilient against various kinds of 
perturbations ||8] |9] [TO] [TT] [12] [T3], temperature fluctuations llT4l and so on. For example one can 
imagine that a strong magnetic field will eventually align all the spins in the direction of the mag- 
netic field and the topological ordered phase transforms to a spin-polarized phase, IfTOl [151 [T6l a phase 
which is easily recognized by local measurements of spins. Or one can imagine that at high enough 
temperature the topological phase transforms to a disordered phase lllTl , again recognizable locally. In 
these transitions a topologically degenerate ground space transforms to other forms of ground states. 
Phase transitions of fliis kind have been studied in|[8l[9l[l0l[IIl[I2[l3l[I2[l6l[Il[l8l[l9l[20l[2ll. 

It is the aim of this paper to study another kind of transition in these models. For concreteness 
we take the Kitaev model and ask how topological order can transform to ferromagnetic order This 
transition is induced by Ising interaction and moreover it signifies a transition between two kinds of 
degenerate ground states. That is in the Kitaev limit the degeneracy comes from topology and in the 
Ising limit, it comes from symmetry. This will certainly add to our knowledge about topological order 
and the way it is either destroyed (i.e. by temperature or by magnetic fields) or changes to other types 
of local order (i.e. by ferromagnetic interactions). 

To this end, we introduce a model in which Ising terms compete with Kitaev interactions, one 
to establish ferromagnetic order and the other to establish topological order We study the model on 
both 2 dimensional torus and the quasi-one dimensional ladder network ||22| which has almost all the 
characteristics of a topological model, i.e. topological degeneracy, robustness and anyonic excitations. 
Studying two different models has the benefit of understanding the role of dimension in this transition. 
To find the possibility of the transition and the transition point if any, we exactly map the problem of 
finding the ground state of the model to a simpler problem. In the ladder case, we map the model to a 
one-dimensional XY model, whose anisotropy is tuned by the ration of Ising to Kitaev couplings. This 
model is exactly solvable by free fermion techniques, its ground state is non-degenerate and smoothly 
varying, except at the extreme points (XX or YY interaction). Therefore in the case of ladder, there is 
no transition at finite Ising coupling. However in the 2D case, we exactly map the problem to the 2D 
Ising model in transverse field. The latter model has been studied using different methods |[23l[24l[25l 
and is known to show a quantum phase transition. We show that the two sides of transition point 
correspond to topological and ferromagnetic order in the Kitaev-Ising model. This provides strong 
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evidence for a transition between these two phases in the original model. 

The structure of the paper is as follows: In section |2]we briefly review some preliminary facts on 
Kitaev model, emphasizing their difference on the torus and on the ladder In section [3] we introduce 
the Kitaev-Ising model and solve it exactly on the ladder in subsection l4. II In subsection l4.2l we map 
the Kitaev-Ising model to 2D Ising model in transverse field and analyze the degeneracy structure 
of the model and interpret it in terms of the original model. In an appendix, by a simple mean field 
analysis, we find the transition point which turns out to be near the actual one obtained by more 
accurate numerical means ll26l |23] l24l |251 . To substantiate the idea of a phase transition between 
topological and non-topological phases, in section |5] we use the above mapping which facilitates an 
estimation of the expectation values of Wilson loops in the two regimes. These estimates indeed 
turns out to be as we expect, that is, the expectation value of a Wilson loop (Wc) behaves as the 
exponential of a quantity which is proportional to the perimeter of the C near the Kitaev point and to 
the area enclosed by C near the Ising point. The paper concludes with a discussion. 



2 A brief account of the Kitaev Model 

In this section we briefly review the Kitaev model 31 in order to set up the notation and use its main 
concepts in the sequel. Consider a lattice whose set of vertices, edges and plaquettes are respectively 
denoted by V, E and P respectively. The number of elements in these sets are respectively denoted 
by |y I, \E\ and \P\ respectively. Spin one-half particles live on the edges of this lattice and hence the 
dimension of the full Hilbert space is given by 21^'. The Kitaev Hamiltonian on this lattice is given 
by 

HK^taev -JJ^^s-KJ^ Bp (1) 

sev p£P 

where 

As := Jl'^*'^' n 

Here i S s means the edges incident on a vertex s and i E dp means the edges on the boundary of 
a plaquette p. The coupling constants, J, K and are taken to be positive. It is easily verified that all 
the vertex and plaquette operators commute with each other and square to the identity operator The 
ground state is thus the common eigenvectors of all the vertex and plaquette operators with eigenvalue 
1, that is I $) is a ground state of the Kitaev model if it satisfies As I <I>) = |<i>), = |$), Vs,p. 



2.1 Kitaev Model on a two-dimensional periodic lattice (a torus) 

For a 2D rectangular lattice of vertices, with periodic boundary conditions on both directions ( a 
torus), we have |y| = A^, ji?! = 2iV, and \P\ ~ N . Hence the dimension of the Hilbert space is 2^^ 
and we also have N vertex operators As and plaquette operators Bp all commuting with each other 
and with the Hamiltonian. However there are two global constraints on the torus, namely 

WAs^I, X{Bp^I, (3) 

.s p 

leading to 2A^ — 2 independent commuting operators and hence a 4-fold degeneracy of the ground 
state. In fact one notes that there are four string operators all commuting with the Hamiltonian, which 
are defined as follows: 

n •= n 
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Figure 1 : (Color Online) The string operators Tl , and , T^?. All of them commute with the pure 
Kitaev Hamiltonian {H{\ = 0)). Only Tl and commute with H{\) for all A. 



tp ■■= n '^^.^ •= n '^^^ 

where Ci and C2 are two homology cycles along the edges of lattice of the torus and Ci and C2 
are two cycles running around the dual lattice (figure 1). Note that these operators, corresponding to 
homology cycles (curves which do not enclose any area) cannot be expressed in terms of vertex and 
plaquette operators. They have the following relations with each other: 

TlTl = -TlTl T^T^ = -T^Tl (6) 

while all other relations are commutative ones. In other words, the operators {T^^Tl) and (T^, T"^) 
form two copies of the Pauli operators a^. and ct^ which act to distinguish the four degenerate ground 
states of the Kitaev model and turn them into each other That is, if we denote the four ground states 
by |$si,s2): '^ii '^2 = 0, 1, then we have 

Tl\^s,.s,) = {-lY'\^s,.s,), T^l*,,,,,) = (-in*,,,,,), (7) 

and 



2.2 The Kitaev Model on the quasi-one dimensional lattice (a ladder) 

Since we will also study the Kitae-Ising Hamiltonian on the quasi-one dimensional systems, it is 
in order to note a few minor differences that the Kitaev model on the ladder has with the 2D case. 
Consider a ladder, as shown in figure (2), with N plaquettes. There are 2A^ vertices and 3A^ edges. 
So the dimension of the Hilbert space is 2"^^. We have periodic boundary condition only along the 
legs. The number of independent operators is equal to N , while the number of independent Ag 
operators is 2N — 1, since A, = /. (There is no such constraint on the operators on the 
ladder). Therefore the total number of independent commuting operators is equal to 3A^ — 1 leading 
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Figure 2: (Color Online) (a) The operators and T.^ on the ladder is the product of ct^ z on the 
lower leg, while is the product of two cr^-'s one on each leg. Both commute with H{X = 0), but 
only commutes with H{X). contains edges from the two legs. No operator containing only 
edges from one leg can substitute it. (b) Electric and magnetic anyons on the ladder 



to a 2-fold degeneracy for the ground state. On the ladder only one pair of operators in (|4] |5)) with 
their properties remain, which are denoted by Tl and in figure (2). In fact the analogue of operator 
is an operator like (Ti ^ sitting on a single rung of the ladder, which no longer commutes with 
the Hamiltonian and the analog of operator is no longer independent from the vertex operators, 
since ~ <Ji,xO'2,x ■ ■ ■ crN,x = A1A2 ■ ■ ■ An, where Ai denote the vertex operators on the upper (or 
lower) leg of the ladder This is in accord with the two-fold degeneracy of the ladder, that is if we 
denote the two ground states of the ladder by I'I's), s = 0, 1, then we have 

T^\<Ps) = ri|<i>,) = |ci>,+i). (9) 

Note that the Kitaev model on a ladder, being a quasi-one dimensional system allows a restricted form 
of topological order That is, concepts like area law for Wilson loops, or topological entanglement 
entropy may not apply to it. However there are still some topological characteristics in the ground 
states. First we have ground state degeneracy which does not come from symmetry. The two states 
being converted to each other by the global string operator T.^. Second we have a finite gap. Also 
the expectation value of any local operator, i.e one which does NOT traverse the the two legs of the 
ladder, is the same on the two ground states \^o) and I'J'i). In fact an operator which distinguishes 
the two ground states, should be one which commutes with the Hamiltonian and at the same time 
anti-commutes with T^. Such an operator is given by which necessarily contains both legs of the 
ladder (It is in this sense that no local operator (one defined on a single leg) can distinguishes the 
two ground states). Finally the system has anyonic exitations of electric and magnetic charges with 
integer charges and abelian statistics. In fact as shown in figure (2), an open string of (Tz operators 
along the edges creates electric anyons at the end points, while an open string of ax operators along 
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the rangs creates magnetic anyons and cycling any electric anyon around a magnetic one creates a 
phase of (-1). Again the specific topology of the ladder reflects itself in the properties of its anyons in 
that, only electric anyons can move around the magnetic anyons. 

Therefore many of the concepts pertaining to topological order are valid also for this quasi-one 
dimensional system. When we speak of topological order on the ladder, we mean this restricted mean- 
ing of the word. On 2D we do not have such a restriction. 



3 The Kitaev-Ising model 

We define the Kitaev-Ising Hamiltonian on any lattice as follows 

H{X) :— Hxitaev + ^Hjsing, (10) 

in which Hxitaev is the usual Kitaev Hamiltonian (1]) and Hjging is the Ising interaction between 
nearest neighbor links 

Hising = — ^ Cri,zO'j,z, (H) 

where {i, j) means nearest-neighbor edges on the lattice. It is important to note that in the presence 
of the Ising interaction, the plaquette operators still commute with the full Hamiltonian, although the 
vertex operators no more do so: 

[Bp,H{X)]^0, Vp, [A,,H{\)]^0. (12) 

Moreover from the four string operators, shown in figure (1), which commutes with the Kitaev Hamil- 
tonian, only two retain this property in the presence of Ising interaction, namely 

[T^,H{X)]=0 [TlH{X)]^0, (13) 

but 

[Ti,i/(A)]^0 [TlHiX)]^0. (14) 

Correspondingly for the ladder, only the operator is defined which commutes with the Hamil- 
tonian. 

In view of the fact that [Bp, H{X)] = 0, Vp, and A and the fact that in both limits (pure Kitaev 
and pure Ising) the ground states have eigenvalue +1 for all Bp's, we conclude that the ground states 
of the Kitaev-Ising model lie in the subspace where Bp ~ 1 for all the plaquettes. Denoting this 
subspace by Vo, 

Vo := {|<^) , Bp\^) = 10)}. (15) 
Therefore the restriction of the Hamiltonian to this subspace, Ho{X) H{X) |vo is given by 

HoiX) = -J^A,-A^t7„(T,,, -if|P|, (16) 

where |P| is the number of plaquettes in the lattice. Therefore the Ising coupling A or more precisely 
the ration j tunes the competition of ferromagnetic order and topological order When this ration is 
zero we have pure Kitaev model and topological order, and when it is very strong, we have ferromag- 
netic order In both cases we have degeneracy, but in one case the degeneracy is due to topology and 
in the other case it is due to symmetry. It is also interesting to note that the degeneracy of the ferro- 
magnetic order is always two-fold, i.e. either all the spins are up or all are down, while the topological 
degeneracy is four-fold for the torus and two-fold for the ladder. In the subsequent sections we will 
understand how this order and the corresponding degeneracy changes as we change the parameter 4- 
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4 Solution of the Kitaev-Ising model 



We showed that the ground states of the Kitaev-lsing model live in the subspace Vg defined in ( fTSl ). 
The restriction of iJ ( A) to this subspace is given by (fTSI l. To further diagonalize Hq (A), we construct a 
suitable basis for the subspace Vq and through this we transform Ho{X) to very simple models which 
have been studied previously. In fact, we will show that for the ladder, Ho{X) turns out to be the 
Hamiltonian of a one-dimensional XY chain, while for the 2D lattice, Hq{\) is the Hamiltonian of an 
Ising model in transverse field. 

The way this basis is constructed is of utmost importance, in fact it should be constructed in such 
a way that all the operators in the Hamiltonian, i.e. the vertex and plaquette and also the Ising terms 
should be represented by nearest neighbor interactions between Pauli operators on virtual spins. Oth- 
erwise, one may come up with an inappropriate reduced Hamiltonian, one which may entail three or 
four-spin interactions or in case of two-body interactions it may entail longer than nearest-neighbor 
interactions. In other words, choosing this basis is a significant step in the process of diagonalization. 
Due to the difference between the topology of the ladder and the torus, we proceed in two different 
ways in construction of this basis. We start with the ladder and then study the case of 2D torus. 

4.1 On the ladder 

Consider the ladder shown in figure (3), where we take for definiteness the number of plaquettes to be 
an even number. We first note that the state 

JV 

|$o):= 11(1 + ^^)1+)^''^' (17) 

i=l 

where |+} is the positive eigenstate of a^, is a ground state of the pure Kitaev model, (one can easily 
check that it satisfies ^^1$) = Bp\^) = for all s and p). Consider the curve Cj on the ladder 
(shown in figure (3)). This is a cycle going around the ladder and in fact it is equivalent to the straight 
curve Ci shown in figure (2) (this equivalence is explained below). Therefore the other ground state 
of the pure Kitaev model on the ladder is nothing but 

:= n ^'.-l*o). (18) 

By equivalence of C[ and Ci we mean that the difference of Jliec' '^i-^ •— IlieCi '^i^^ ^ 

product of Bi operators which has no effect on |$o). In fact we can simply straighten a □ or a U by 
multiplying with the B inside them. We now construct the following set of un-normalized states 

2N 

|f ) := f2, fa, • • • f2jv) = n <:J*o), n = 0, 1. (19) 

i<£C[ 

Clearly these states satisfy i?p|f ) = |f ) for all p. We also note that 

|{r~,;=0}) = |$o) , |{f,: = l})-|$i). 

Moreover, they are orthogonal. For the proof of orthogonality, the basic idea to use, is that |$o) can 
be viewed simply as a linear combination of closed loops of spin — particles in a background of all 
spin + particles. Now if r 7^ r', then it is easy to see that (f |f') = ($o| riiGC '^i'z \^^) '^^e 
product of two states where one (I'I'o)) has all closed loops and the other has open strings of spin 
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Figure 3: (Color Online) The curve used for generating an orthonormal basis in the sector Vq for the 
ladder network. 



down states, the product of which is zero. Finally their number is 2 which is equal to dimension of 
Vq. Thus with proper normalization, they form an orthonormal basis for Vq. Why we have chosen 
this particular curve and this particular form for expressing the states of this sector? The answer lies 
in the nice form (i.e. nearest-neighbor two-body interaction) of the reduced Hamiltonian Hq{X). If we 
choose the curve as a simple straight form like Ci, then the states |r) do not span the whole subspace 
Vo. 

To find Hq{\) in this new basis, we should determine the action of operators Ag and also the 
Ising terms on these basis states. Due to the zigzag shape of the path and the appearance of the 
corresponding operators in the definition of |f ), we find that any vertex operator like Ai in figure (3) 
when acting on the state (T% passes through all the operators except (Tza-i and ^ with which it 
anticommutes. Thus the passage of Ai through the whole chain of operators produces only the factor 
( — hence the following effective operation on the basis states: 

A,|f) = Z,_iZ,|f), (20) 
where Zi is the notation of Pauli operator in this subspace. 

Notation: Original Qubit states on the edges of the lattice are denoted without a~. Thus |0) and 
|1) denote the computational basis states on the edges, |0) = \z, +) and |1) = |z, — ). Pauli operators 
on these qubits are denoted by ax and cr^. The qubit states in iT% are always denoted by a~and the 
corresponding Pauli operators on these qubits are denoted by capital letters, X and Z. 

We now come to the Ising terms. Consider a group of Ising terms in one plaquette, the one shaded 
in figure (3). This can be written as 

Ci := ai^zO'i+l.z + O'i^zO'a.z + 0'i+i^zO'i+2,z + 0'a,z(Ti+2,z ■ (21) 

To express the action of Ci in the basis ( fT9] l, we note that this can be rewritten as 

C\ := {(Ti,zO'i+i,z + (^i+i,zO'i+2,z) (1 + Bj), (22) 

where Bi is the operator corresponding to the same (shaded) plaquette. The operator (1 + Bi) gives 
a factor of 2 when acting on the state |f ) and the remaining <Tz.i operators only flip the corresponding 
bit labels f^. Hence the following effective action on the state: 

C,|f) = 2 iX,X,+i + X,+iX,+2) |f), (23) 
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Figure 4: (Color Online) The state Oq, defined in equation (l25T l and the following paragraph. 

where again Xi is used to denote the first Pauli operator on the subspace Vq- Putting everything 
together, we arrive at the following effective Hamiltonian on this subspace: 

HoiX) = -JY,Z,Z,+i - 2\Y,X^X,+l - KN. (24) 

i i 

This is an XY Hamiltonian in the absence of external magnetic field which has been studied ex- 
tensively in the literature ||27i . Its exact solution is provided by turning it into a free fermion model 
by Jordan-Wigner and Bogoluibov transformations. Its ground state is non-degenerate except at the 
extreme points A = or A — > oo. The ground state and the correlation functions show no non- 
analytical behaviour and no quantum phase transition occurs for finite A. The only thing which hap- 
pens is that the two-fold degeneracy breaks for any value of A except at the extreme points A = 0, 
(Pure Kitaev) or A — > oo (Pure Ising). The end conclusion is that a transition from topological to 
ferromagnetic order does not occur for finite A in quasi-one dimensional systems. 

Before leaving the subject of ladders, it is instructive to have a final look at the ground states of 
Ho{X) at the two extreme points. In these two limits, the ground state(s) of ( l24l i should have simple 
product form (in terms of the labels fi). Let us see if these are really what we expect for the Kitaev- 
Ising model. In the limit A = 0, equation (l24l i says that the virtual spins should all align either in 
the positive or negative z direction, hence there are two degenerate ground states given by |0) and 
|1) respectively. As explained at the beginning of this subsection, these two states are clearly the two 
Kitaev states |$o) and |<I>i) on the ladder The other limit, however, is more tricky to show. In the 
limit A — > oo all the virtual spins should align either in the positive x or negative x direction. We 
should show that this means that the actual spins on the edges of the ladder all align in the positive or 
negative z direction. To see this consider the state |+, +, ■ • ■ +). In view of the definition ( fT9] l and the 
structure of (fTTI l. and the fact that |+) cx |0) + this corresponds to 

N 

I+, +,■■•+)= l[il + a,^.)\'i'o)^l[{l + Bi)\no), (25) 

where jfio) := Iliec' ~^ <^i,z)\+)^^^ is a state on the ladder which we depict in figure (4). Here 
we have used the property (1 + (Jz)\+) = |0). When the operators 1 + Bi act on \^lo), they turn the 
remaining + states into and hence turn it into a ground state of the pure Ising model. 

For the other state |— ) a similar reasoning works in view of (1 — (Tr)|+) = |1) where |1) is spin 
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Figure 5: (Color Online) A portion of the lattice with their Kitaev and Ising interactions. The Ising 
interactions are between nearest-neighbor links. An Ising interaction like az.zCJ^.z (shown with the 
dash line) commutes with all the vertex and plaquette operators except with Aa and Aa' , to which it 
anti-commutes. Similarly an Ising term like a2.zO'3,z commutes with all but Ai, and Ai,r. 

down in the z-direction. 



4.2 On the two dimensional lattice 

We now turn to the square lattice and show that Ho{X), describing the interactions of virtual spins, 
is in fact the Hamiltonian of a 2D Ising model in transverse magnetic field. This model is known 
to undergo a transition from ferromagnetic order to spin-polarized ordered phase. These phases, are 
shown to correspond respectively to topological and ferromagnetic ordered phases for the actual spins 
on the edges of the lattice. 

To show this equivalence, we follow steps similar to the ones in previous section, however to 
represent the Hamiltonian in a simple form, we should choose an entirely different basis for the sub- 
space Vq. Let the lattice have N plaquettes. Then the number of edges will be 2N and the number 
of vertices will be N. In a concise notation we have \P\ = N, \E\ = 2N, and \V\ = N. The 
dimension of the full Hilbert space is thus 2^^. Vo is the common eigenspace of all Bp operators 
with eigenvalue +1. Since the number of independent plaquette operators is iV — 1, this means that 
dim{Vo) = 2!^-i = 2^+^. Furthermore this subspace is decomposed to four different disconnected 
subspaces according to the eigenvalues of the global string operators and T^. Let us denote this 
decomposition by 

Vo = Vo++ ® Vo+- © Vo"+ e Vq— . (26) 
Each subspace is 2^^^ dimensional. Consider now the states 

If) n^riO)^""' n^Q,l. (27) 
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Obviously these states satisfy Bp\r) = |f ). Moreover, when an Ai acts on these states, it increases (by 
mod 2) the label r^, hence the action of each Ai on these states is represented by the bit-flip Palui op- 
erator Xi. In view of the constraint IljeB ^ 1' have the equality |f ) = |f ), where = ?'i + 1, 
mod 2 Vi. The subspace Vq"^ is therefore span of the equivalence class of states [|f)] = (|f), |f)). 
The other subspaces can be constructed similarly, i.e. Vq^ = span{[T^\r)] = (T^|f),T^|f))}. 

When A — > oo, (pure Ising model), the two ground states of the pure Ising model are clearly in 
the subspace Vq^+, hence by the fact that [T^,H{X)] = [T^,H{X)] = and by continuity we find 
that the ground states of H{X) live in the subspace Vq'^. Hereafter we will focus on this subspace. 

To proceed we also assume that the lattice is bi-partite, i.e. V = Va^ Vb, where the vertices in 
Va are denoted by black circles in figure (5) and those of Vb are denoted by white circles. Note that 
this puts a condition of even number of vertices in both directions. We need to find the action of the 
operators Ai and the Ising terms on the states (|27] |. It is obvious that the action of a vertex operator 
like Ai on the state dZTl i is to simply flip the bit r^, therefore Ai acts on this subspace as Xi. Next we 
come to the Ising interactions. Consider the shaded plaquette in figure (5). The Ising interactions are 
given by 

0'l,zO'2,z + 0'2,zO'3,2 + a3^z(^i,z + Cri,z<^l,z = {(^l,z<^2,z + <^i,z<y 4,,z){'^ + B) , (28) 

where B is the plaquette operator containing the links 1, 2, 3, and 4. We now use the fact that an 
Ising interaction like 0'i.zcr2,z commutes with all the vertex operators and anit-commute with Aa and 
Aa'- This means that when cri,zCT2,2 acts on the state |f) it simply produces a factor that 
is this operator acts on the subspace Vq'"^ as ZaZa'- Similarly the Ising term (Jx.zcri^z commutes with 
all the vertex operators and anti-commutes with A^ and At,' and with the same reasoning the action 
of this operator on V^^ is equivalent to Z^Ziji. Therefore the Ising terms couple the nearest neighbor 
vertices of the sublattice Va and Vb separately. Putting everything together we find the following 
effective Hamiltonian: 

H^Ha + Hb- KN, (29) 

where KN comes from the action of Bp on Vq and Ha and Hb are each a 2D Ising model in 
transverse field on sublattice A and B respectively; 

is A l^i^j)tzA 

and 

Hb = -JY.X,-2\ ^^^r (31) 

Here (z, j) means nearest-neighbor vertices on the corresponding sublattice. Note that the factor 
of 2 in front of A comes from the factor (1 + B) in ( |28] |. 

In this way the Kitaev-Ising Hamiltonian turns into the rather well-studied 2D Ising model in 
transverse field. The ferromagnetic order is controlled by the Ising coupling which tries to align all 
the virtual spins in the +z or — z direction. The transverse magnetic field controlled by J competes 
with the Ising interaction and destroys the order if J passes a critical value Jc- Density matrix renor- 
malization group ll24l gives a value of the critical magnetic field as Jc ~ 6A. A simple mean field 
analysis (provided in the appendix) gives the value Jc = 8A. When J > Jc, the virtual spins try to 
align in the +x direction. It is important to note that in the limit J = (or A — > cxd), the ground state 
of the virtual spin system is doubly degenerate, while in the limit J — > cx) (or A = 0) the ground state 
is unique and non-degenerate. What is interesting is that these two phases of virtual spins correspond 
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to the topological and ferromagnetic phases of the actual spins on the edges of the lattice. 

To see this correspondence, consider one of the Hamiltonians, say Ha- The J term tends to align 
all the virtual spins in the +x direction, while the A term tends to align them in the positive or nega- 
tive z direction. In the limit A — > we have a unique ground state |+, + , • • • +), while in the limit 
A — > oo, there are two ground states |0, 0, • • • 0} and |1, 1, • • • 1). The same thing happens in sub- 
lattice B. Therefore in the limit A — > 0, there is a unique state, denoted by (-1-^, +b), which is the 
topological ground state of the Kitaev state |$^^), while in the Umit A — > oo, there are two ground 
states (0^, Ob) = (1^, 1b) and (Oa, Is) = (1a, Os) which are the ferromagnetically ordered states, 
where all virtual spins are either up or down in the z direction. Let us show this in a more explicit way. 
Consider the state denoted by {+a,+b)- In view of the notation (|27]|. and the fact that \+) oc |0) + |1) 
the state of actual spins corresponding to this state is given by l^"*"^) := ni(Ev(^ which 
is a common eigenstate of all the Ai and Bi operators with eigenvalue 1, hence a ground state of the 
pure Kitaev model. 

In the other limit, the ground state {Oa,Ob) denotes the state (IZTT i. where none of the Ai's act on 
the state lO)*^'^', hence this is nothing but a uniformly ordered ferromagnetic state in which all the 
spins are up in the z direction, i.e. lO)*^'^'. We remind the readerthat (Oa, Ob) = (Ia, Ib), due to the 
constraint Jlie v — 1- other words, if we act on the state jO)®'^' by all the vertex operators Ai, 
nothing happens since the flipping actions of Ai operators on the sublattice A are neutralized by those 
on sublattice B. To flip all the spins, one needs to apply the vertex operators A, on only one sublattice, 
hence the states (Oa, Ib) = (1a, Ob) correspond to the ferromagnetically ordered state jl)®'^'. 



5 Topological characteristics; estimates of Wilson loops 

In order to justify the transition from topological to ferromagnetic order, we can estimate the value of 
a Wilsonian loop, 

(Wc) ■■={l['y^,.), (32) 

where the expectation value is calculated in the ground state and C is a closed curve on the dual lattice 
enclosing an area 5, i.e. C = dS. Let us denote the perimeter of dS by | dS \ and the area of 5* by 151. 
Then it is known that in the topological phase, the expectation value of this Wilson loop behaves as 
while in the non-topological phase it behaves as e'^'I'^l, where (3 and 7 are two constants 1281 . 
The mapping of the Kitaev-Ising model to the 2D ITF model allows to obtain estimates of the Wilson 
loop in the two regimes perturbatively. To proceed we first note that in view of the definition of the 
vertex operators As in (O, the operator Wc can be written as 

Wc = Y[As = l[Xs, (33) 

where in the last equality we have used the equivalence of As with Xs on virtual spins (see the 
paragraph after Eq. l27t . where the Hamitonian becomes a simple ITF Hamiltonian as in ( l30l ). we have 
to calculate the following expectation 

{Wc) = ^ (34) 

where \^) is the ground state of the 2D ITF model with the Hamiltonian given in ( |30] l. Note that in 
view of the decoupling between the two sublattices, we only consider one sublattice. Consider now 
the two limits, near-Kitaev and near-Ising separately. 
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5.1 Close to the Kitaev limit 



In this limit, where 7 := j << 1, we can take the Ising term in ( [30] l as a perturbation to the magnetic 
field and approximate the ground state j^E') as a series 

I*) = 1*0) + 71*2) + 7'(l*2) + 1*4)) + • • • . (35) 

Here |^'o) = |+)®^ is the ground state in the limit A = and |^'2) denotes the linear superposition of 
all states in which two adjacent spins have been flipped by the ZiZj terms. Note since we are doing an 
estimate and also we do not assume these states to be normalized, all numerical factors coming from 
perturbation expansion like energy differences and so on are absorbed in the definition of these states. 
Similarly, \^a) is the linear superposition of all states in which 4 (two pairs of nearest-neighbor) spins 
have been flipped due to {ZiZj){ZkZi) terms and \ '^2) is the Unear superposition of all states in which 
only two non-adjacent spins have been flipped by {ZiZj){Zj Zk) term and so on. We then have 

(^-l^-) = (^-ol^o) + 7'(*2|«'2> + 0(73). (36) 
We now note that |\l/2) can be broken up into three kinds of states, i.e. 

|*2) = |vI/2)^+|*2L + |*2)os, (37) 

where these states are described in figure (6). In view of this figure and the fact that Xi\±) = ±|±), 
we then have 

W^C|*2) = 1*2)^+1*2)^ - l*2)as- (38) 

Combining dJTl l and ( l38T l. and keeping all the terms up to order 7^ 

2H|ff)=l-2.^..(*.|*.)„.+Oh») (39) 
and since (^'2|'I'2)as ^ \dS\, we find that close to the Kitaev limit, we have 



la^l. (40) 



Therefore as expected close to the Kitaev limit, the expectation value of the Wilson loop behaves as 
the exponential of the perimeter of the loop, which is characteristic of the topological phase. 



5.2 Close to the Ising limit 

We now consider the Ising limit where j := ^ < 1- Right at the Ising point, consider one of the 
degenerate ground states, say one in which all the spins are in the +z direction, or using the quantum 
computation terminology, all the spins are in the state 0. The magnetic field in ( [30l ) perturbs this 
uniform ground state by flipping spins one by one. These spins can be inside and or outside the loop 
dS. Denote the lattice points inside the loop by S and the lattice points outside it by 5*. Let us \(pk)\Xk) 
denote the product state in which \<f>k) is the state pertaining to S, and is the uniform superposition 
of all basis states in which exactly k spins have been flipped to |1) (or | — z)) and |xfe) is the state 
pertaining to S in which any number of spins have been flipped. Therefore the state \4>k)xh) is a state 
in which at least k spins have been flipped. Then we can write the perturbative ground state as 

fe=0 
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(a) 



(b) 



(c) 



Figure 6: (Color Online) Different contributions to the state |^'2)- All the spins are in the state |+), 
except the two spins at the end of the bold link which have been flipped by the ZZ interaction and 
are in the state | — ). Figures (a), (b) and (c) exemplify contributions to the states \^2)s, 1^2)5 and 
1*2)95, respectively. 

Note that \xk) is a state which is normalized to 0(1). The reason comes from perturbation theory, 
that is \xk) is the superposition of states in S in which 0, 1, or more spins have been flipped. Note 
also that 




Since the operator Wc = Hies -^"^ '^^P^ ^ spins inside S, we have 



\s\ 




(42) 



Therefore we find 



($1$) 




(43) 



On the other hand we find from (|42] | that 



(0A;|0fe> (XfclXfe) 




(44) 



Dividing (l44l i by ( |43] | we find 




(45) 
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Therefore we have shown that close to the Kitaev and the Ising points, the Wilson loop behaves as 
expected, that is, its logarithm is proportional to the perimeter of the loop in the topological phase and 
proportional to the area in the ferromagnetic phase. All this has been made possible by mapping the 
system to the 2D Ising model in transverse field. 

6 Discussion 

We have introduced the Kitaev-Ising model, equation ( fTOl ) as a model for studying the transition be- 
tween topological order and ferromagnetic order in a lattice system. In particular we have shown that 
on the quasi-one dimensional system of the ladder (with periodic boundary condition), there is no 
quantum transition between these two kinds of order at finite A, while in two dimensions a transition 
occurs at finite A. This is reminiscent of what we have for thermal phase transitions based on symme- 
try breaking of discrete symmetries. 

In the quasi-one dimensional case, we have exactly mapped the problem to the problem of find- 
ing the ground state of an XY chain in zero magnetic field for which exact solution by free fermion 
techniques is available. On a two dimensional lattice on the other hand, we have mapped the ground 
sector of the Kitaev-Ising Hamiltonian to two copies of Ising models in transverse magnetic fields, 
each defined on one sublattice, the latter model known to show sharp transition for finite A. Although 
we have not attempted a numerical study of the model near the transition point, the equivalence with 
the 2D ITF model combined with the analysis of the degenerate structure of the ground states and 
their global properties in the two limits show that such a transition does occur for some finite A. We 
have also estimated the Wilson loops and have shown that close to the Kitaev and Ising points, the 
logarithm of the expectation value of a Wilson loop is proportional to the perimeter of the loop in the 
topological phase and to the area enclosed by the loop in the ferromagnetic phase. It is also worth 
noticing an intriguing difference between the characteristic of the two different phases. In the topo- 
logical phase, the four ground states are distinguished by loop operators and and are mapped 
to each other again by loop operators and T^. In the ferromagnetic phase on the other hand, the 
two degenerate ground states are distinguished by a local operator cr^, while the two ground states 
are mapped to each other by a global operator HieB '^i-x^ encompassing the whole lattice. Therefore 
during the transition, the distinguishing loop operators shrink to points, while the transforming loop 
operators expand to the whole lattice. 

This study can be extended in a few directions. First, one can use numerical techniques to deter- 
mine the ground state and its properties as a function of the Ising coupling. Second it is desirable to 
generalize the analysis of this paper to the cases where the 2D lattice is not bipartite or the number 
of plaquttes in the ladder is not even (the simpUfying assumptions made here) and to see if it leads to 
significantly different results. 
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8 Appendix 

In this appendix we briefly do a mean field analysis of the 2D Ising model in transverse field. Such an 
analysis reveals only a very qualitative feature of the transition. Using a product trial wave function 
1^) = 10)^^ for ths 2D ITF model dSOl l, one needs to minimize the energy 

:=-J(,/)|X|0)~4A(0|Z|0)2. 

Taking |0) = cos ||0) + sin |e*'^|l), leads to the following expression 

e(6l,0) = -Jsin 6*008 </>-4Acos2 61. (46) 

Minimizing this energy, one obtains that the nature of the mean field ground state changes at a critical 
value 7c = 1 ( 7 := ^), that is the state which minimize the mean field energy is 

I COS -^|0) ± sm ^|1) t < 

where 

sm9^ = ^. (48) 

From this mean field analysis we find the following expectation values {X) = {As) and (Z) = (ctz), 
also shown in figure (7). 

{X)^{As)^l^} ]<J^ (49) 



8A 



and 



A < 



J 
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Figure 7: (Color Online) Mean field expectation values of {X) and (Z) on the 2D Ising model in 
transverse field, which are respectively equal to the {As) and (az) on the original lattice. Note that 
after the transition point, the ground state becomes degenerate and the two ground states can be dis- 
tinguished by a local order parameter (az), hence two branches. The transition point predicted by 
mean-field is perturbation theory ||26) and renormalization method ll24l give a value « jr. 
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